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STUDY REPORT 


THE APPLICATION OF SYSTEMS ANALYSIS AND 
MATHEMATICAL MODELS TO THE STUDY OF 
ERYTHROPOIESIS DURING SPACE FIIGHT 




1, 0 INTRODUCTION 


It has been demonstrated that individuals returning from space missions - 
from 4 to 85 days - have exhibited a significant loss of their preflight red 
blood cell mass. Some of the various factors that have been implicated in 
attempting to explain this loss include weightlessness, immobility of crew, 
high atmospheric oxygen partial pressure and lack of atmospheric nitrogen 
(1, 2, 3). A universally agreed upon explanation for this phenomena is still 
lacking. 

This study is designed to be a preliminaiy systems analysis of the physio- 
logical systems that control red cell production. One of the purposes of 
this report is to review the physiological mechanisms that are presently 
believed to control erythropoiesis. A second objective is to evaluate existing 
mathematical models of erythropoiesis regulation. With this as a basis it 
was hoped to formulate a new model - either a composite of existing models 
or one embodying improved features - that would form the starting point for 
simulating the erythropoietic response to weightlessness. Finally, an assess- 
ment would be made regardir^ the role that systems analysis can play in the 
evaluation of the Sliylab hematologic program. 

Systems analysis has been shown to be a highly useful technique in helping to 
understand complex physiological processes, especially feedback control 
systems. This method usually consists of several phases of development: 

1) the formulation of mathematical models of these processes, 2) a study 
of the behavior of these models using computer simulation techniques, and 
3) the comparison of the behavior of the idealized model with the real system. 
The development of the model should be based, wherever possible, on accepted 
experimental evidence and physiological concepts. It is usually the case, how- 
ever, that a complete understanding of the system does not exist. In those 
situations, systems analysis can be vtsed to test various hypotheses quickly 
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and efficiently. TMs method finds its greatest utility when combined with a 
positive program of physiological research aimed at clarifying the hypothesized 
events. 

2. 0 PHYSIOLOGY OF ERYTHROPOIESIS 


Factors that Control Erythropoiesis 

Our present understanding of erythropoiesis, although far from complete, is 
sufficiently evolved to define a sequence of mechanisms that describe the 
regulation of red cell production (5, 6, 7, 8). It can be seen from Figure 1 
that such a sequence forms a closed feedback loop in which red cell produc- 
tion is adjusted according to the balance of tissue ojjygen supply and demand. 
There is general agreement that the direct stimulus for red cell production 
in the bone marrow is an erythropoietic stimulating factor known as ESF or 
erythropoietin. Release of erythropoietin, thought to take place primarily 
from some "oxygen detector" tissue in the kidney, is regulated by the balance 
between oxygen supply and oxygen demand. Decreasing the oxygen supply in 
relation to the demand results in increasing the rate of erythropoietin release 
and vice versa. Erythropoietin is released into the general circulation and its 
final plasma concentration is determined by its rate of release, volume of 
distribution and the rate at which it is metabolized. The target organ of 
erythropoietin is the bone marrow where it is believed the ESF plasma concen- 
tration determines the production rate and release of red cells. Although the 
mechanism of action has yet to be clarified, er 3 dhropoietin is thought to pro- 
mote the mitotic activity of er 5 dhroid stem cells as well as stimulating hemo- 
globin production. Etythrocytes, once released and distributed into the cirr- 
culation, become oxygenated in the lungs. The blood concentration of oxyhem- 
oglobin and flow rates are the major determinants of the degree of tissue 
oxygenation which determines the rate of cell production. This oversimplified 
view accounts for the most important elements of the eiythropoietic system. 
Some of the details of this system will be discussed below. 
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FIGURE I 


HYPOTHESIS FOR CONTROL OF ERYTHROPOIESIS 





Physiology of Tissue Oxygenation 

The amount of oxygen deliyered to a tissue depends on a number of factors; 
oxygen tension of inspired air, pulmonary function, hemoglobin concentra- 
tion, affinity of hemoglobin for ojg'^gen, speed of dissociation of ojqrgen from 
hemoglobin, cardiac output, and the vascular distribution of circulation blood 
among the individual tissues. From this amount of oxygen, a certain fraction 
is extracted and consumed by the tissue. Oxygen enters the cellular spaces 
by diffusion along an oxygen tension gradient between the capillary and the 
cell, A change in any of the factors listed above or in the cellular metabolic 
rate will result in a change in tissue oxygenation. 

Thus, there are many factors in addition to hemoglobin concentration that 
could cause a temporary imbalance between tissue oxygen supply and demand. 

Any small change in this balance is usually compensated for by rapid adjust- 
ments in the circulatory and respiratory systems. For example, changes in 
the number of active tissue capillaries and blood flow adjustments are the 
most common compensatory responses to minor local fluctuations in the degree 
of tissue oxygenation. Under these conditions the rate of red cell production 
and the circulating red cell mass are remarkably constant. However^ in more 
extreme and chronic situations such as hemorrhage, altitude hj^ioxia, or pul- 
monary disorders, an increase in red cell production appears the preferred 
way by which the body compensates for the resulting insufficient oxygen trans- 
port to the tissues. In addition, conditions that lead to extreme levels of 
tissue hyperoxia such as breathing from an hyperoxic atmosphere or trans- 
fusion induced polycythemia or when the tissue demand for oxygen is decreased 
as in hypothyroidism, hypophysectomy, and starvation there is a significant 
decrease in erythropoiesis. 

Available evidence strongly supports the presumption that changes in cell production 
and circulating red cell mass are triggered by changes in the tissue oijQ^gen tension. 
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But, as we have discussed, the circulatory system is also capable of respond- 
ing to a change in tissue oxygen tension. Hence, it is important to understand 
the delicate balance that may exist between the responses of the circulatory 
and erythropoietic systems. To say that the circulatory system responds to 
minor fluctuations in tissue o^grgen tension and that the erythropoietic system 
responds to more extreme and prolonged changes is no longer adequate, 
especially if one is attempting to construct an integrated model describing 
overall circulatory function. It is quite possible, although not yet proven, 
that the erythropoietic system is temporarily activated by even minor imbal- 
ances between ojQrgen supply and demand. On the other hand, the character- 
istics of the still unknown tissue oxygen sensing sites may be such that a 
sufficient time delay exists before erythropoietin or its precursors are released. 
This could allow time for circulatory changes to take place that may fully com- 
pensate an imbalanced condition without stimulating the slowly responding 
erythropoietic system. A more complete understanding of this problem may 
have to await the identification and characterization of the special sites that 
are believed capable of responding to changes in tissue oxygenation, 

A decrease in oxygen-hemoglobin affinity resulting in improved oxygen release 
can be considered comparable to increases in blood flow or hemoglobin con- 
centration in promoting the supply of oxygen to tissues. Changes leading to 
this shift to the right of the oxyhemoglobin saturation curve have been observed 
in subjects exposed to hypoxic conditions at altitude, exercising individuals and 
anemic patients. A constituent of the red cells has recently been implicated 
in this phenomena when it was discovered that the concentration of 2,3- 
diphosphoglyceric acid (2, 3-DPG) was abnormally high in these situations (8). 
The dissociation of oxygen appears to be facilitated by 2, 3-DPG because it 
preferentially and reversibly binds deoxyhemoglobin. Conversely, conditions 
of relative hyperoxia, such as moving adapted individuals from altitude to sea 
level, has been observed to be accompanied by decreases in the levels of 
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2, 3-DPG. It is understood that the blood analysis of the Skylab crewmembers 
includes an analysis of 2, 3-DPG. This data might be useful in estimating the 
relative shifts of the oxyhemoglobin curves both during flight and upon recov- 
ery, and incorporating this information into an integrated model of the erythro- 
poietic system. 

Renal Release of Erythropoietin 

The majority of experimental evidence reviewed in recent years (7, 8, 9, 10) 
supports the hypothesis that ojq/^gen supply (as measured by tissue PO^) to a 
renal ’’detector” site in relation to the oxygen demand of that region is the 
primary stimulus for erythrocyte production via the release of a renal erythro- 
poietic factor (ESP). However, the renal site of origin or activation of the ESF 
has not yet been found. It has been suggested (7, 10) that the kidney tissue has 
developed a specialization that enables it to function as a sensitive 
ojqrgen chemoreceptor. While most organs in the body exhibit a blood distri- 
bution pattern that results in an even release of blood oxygen, the kidney has 
a peculiar microcirculation: a) the sharp angles between arterioles and 
arteries leads to plasma-skimming where the blood in some of the smaller 
vessels have a reduced hematocrit; b) efferent arterioles from the cortex dip 
deep into the medulla and then return as capillaries, resulting in a steep grad- 
ient of tissue oxygen tension along this axis. Thus, a relative hypoxia is 
present which may act as a continual stimulus for daily erythropoietin produc- 
tion. Any small change in blood oxygen delivery should have a marked effect 
on this area of low tissue PO^. 

Although a large nmnber of experiments have shown that changes in red cell 
production are accompanied by corresponding changes in plasma ESF, this 
has not always been observed. In order to explain some of these "unexpected" 
results it has been speculated that there are agents formed in the plasma that 
act to accelerate or inhibit er 5 d:hropoiesis independently of ESF. However, 
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others (7) seem to feel that many of the apparent discrepancies appear to 
result from the low level of resolution of various bioassays for erythropoietin 
and the indirect means of estimating total oxygen supply. They do not believe 
that existing data call for the postulation of additional chemical factors. Thus, 
it would seem that the theory of oxygen supply and demand is the best explana- 
tion we have for the existing data. 

Even if an additional system of erythropoietic regulation besides ESF does not 
exist, there is a growing body of evidence to suggest that ESF is not formed in 
the kidney, but that the kidney produces a renal erythropoietic factor (REF) 
which can be activated to ESF by unknown constituents in normal serum. If 
this precursor to ESF does exist, the dynamics of the erythropoietic system 
could follow the laws of second order, rather than first order, kinetics, 
mathematical simulation of this system might be improved by including this 
effect. 

Renal blood flow could markedly affect the amoimt of oxygen delivered to the 
special sensing sites. Moreover, in contrast with all other body organs and 
tissues, the oxygen demand or consumption in the kidney is more or les s 
proportional to blood flow rather than remaining constant (11). '^’hus, while 
increasing blood flow may increase oxygen supply there is a proportional 
increase in oxygen demand. In addition, renal blood flow is known to be under 
a high degree of autoregulatory control., These considerations suggest that 
the balance between tissue oxygen supply and demand is not markedly influenced 
by blood flow in the kidney. Whether the renal oxygen detector sites exhibit 
these same characteristics that may exist for the kidney as a whole is not known. 
Several experimental studies have shown that reduction in renal artery blood 
flow resulted in significant increases in plasma levels of erythropoietin (e. g. 12). 
However, in summarizing the total body of evidence for the effects of renal 
blood flow on er 5 d.hropoiesis, several reviewers have appeared to reach 
opposite conclusions (7, 10). 
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3.0 ' COMPARISON OF MODELS FOR THE CONTROL OF ERYTHROPOIESIS 

A review of the published literature reveals only three models useful for math- 
ematical simulation of the er 5 d;hropoietic control system (5, 6, 13, 14). The follow- 
ing section is a review of these models discussing their similarities and differ- 
ences. Comparison will he facilitated by breaking each model into component 
parts as suggested in Figure 1. Figures 2, 3, and 4 are flow charts of the models 
of Mylrea and Abbrecht, Hoc%son, and Guyton*. An explanation of symbol defini- 
tions used in Figure 4 appears in Table 1. Figure 5 is a summation of the dis- 
cussion to follow. 

Oxygenation of Hemoglobin {2/1&2; 3/4; 4/3 )** 

The degree of oxyhemoglobin saturation in arterial blood is simply an input term 
in the Hodgson model. Mylrea and Abbrecht account for changes in blood PO^ in 
response to acute and acclimatized hypoxia using an empirical relationship and 
then compute hemoglobin saturation from an experimental oxyhemoglobin satura- 
tion curve. The percent hemoglobin saturation is computed by Guyton from a 
pulmonary dynamics section of his model and appears as an input term in the 
red cell production block of Figure 4. 

Oxygen Delivery to Kidney (3/2; 4/1) 

Mylrea and Abbrecht consider blood flow to the kidney to be constant and do not 
include a blood flow term in their model. Hodgson utilizes renal blood- flow as , 
an input variable, but also states that "it is not kidney blood flow that is the sig- 
nificant factor, but the flow to the (renal) tissue in which the (renal tissue partial 
pressure) detector is located". Guyton considered the site of erythropoietic stim- 
ulation to be located in "non-muscle, non-renal tissue". Consequently, the blood 
flow term that is used in his model (computed from a circulatory model) is non- 
renal blood flow which appears as an input variable in Figure 4. 

* Figure 4 is a portion of the complete Guyton model (13jl4) concerned only 
with the control of red cel 1 production. This model will become the basis 
for a proposed model for future study. There are some minor differences 
between Guyton'S' section for red cell control and the representation in 
Figure 4. These will be explained more fully in the text. 

** This notation refers the reader to the appropriate flowchart and block number. 
In this example, oxygenation of hemoglobin is described in Figure 2, Blocks 
1 and 2; Figure 3, Block 4; and Figure 4, Block 3. 
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FIGURE 2 


MODEL FOR CONTROL OF ERYTHROPOIESIS 


Erythropoietin 

Kidney distribution volume Bone morrow 



hemoglobin concentration, Vi = viscosity factor, HbO = eifective oxyhemoglobin 
concentration, R = rate of erythropoietin release, E = plasma erythropoietin con- 
centration, Eq = normal plasma erythropoietin concentration, = distribution volume 
for erythropoietin, P = rate of hemoglobin production, MT -- erythrocyte maturation 
time, L — rate of hemoglobin loss, TH = total circulating hemoglobin, [Hb] = blood 
hemoglobin concentration, Fb = blood volume, Fp — plasma Volume, Fpo == normal 
plasma volume, Fpr = steady state hypoxic plasma volume, MCV — mean corpuscular 
volume, MCH = mean corpuscular hemoglobin, k — constant. 

(MYLREA& ABBRECHT, 1971) 
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FIGURE 3 


MODEL FOR THE CONTROL OF ERYTHROPOIESIS 
(HODGSON, 1970) 



CHb 02 — oxygen capacity of Hb, EP = total erythropoietin, ( EP) = plasma erythropoietin 
concentration, FUO 2 = fraction of maximum oxygen flow delivered to renal tissues, (Hb) = 
hemoglobin, concentration, MCHb = mean corpuscular hemoglobin, = oxygen tension 

of renal tissue, = blood volume, Q|^(t) = blood flow rate, = red cell volume, 

R*^C = red cells produced daily, RC = red cell count, = oxygen saturation of art. 

blood, = oxygen saturation of ven. blood, Tep = mean life span of erythropoietin, 

7 = mean life span of red cells, = maximum oxygen flow rate in arterial blood, 

VO 2 == tissue oxygen consumption rate. 
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FIGURE 4 


GE MODEL FOR CONTROL OF ERYTHROPOIESIS 
(DERIVED FROM GUYTON) 


B 

L 
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0 

D 



11 






TABLE I 


SYMBOLS AND DEFINITIONS OF TERMS USED IN FIGURE 4 


GUYTON 

this 



NOTATION 

REPORT 

DEFINITION 

UNITS 

AOM 

AOM 

Autonomic stimulatory effect on tissue oxygen utilization 

ND 

*BFN 

BF 

Renal blood flow 

i/min 

VB 

BV 

Blood volume 

liters(i) 

*02M 

BM02 

Basal oxygen’ utilization in renal tissue 

ml/ min 

** 

CHbOg 

Oxygen carrying capacity of hemoglobin 

ml/gm 

OVA 

002 ^ 
art 

Volume concentration of oxygen in arterial blood 

ml/l 

P02 

DF 

Drive factor causing red cell production 

mm Hg 

- 

F02, . 
tis 

Fraction of maximum oxygen flow delivered to tissues 

ND 

HM 

HOT 

Hematocrit 

ml/lOO ml 

POY 


Sensitivity of red cell production 

i/mln-mm Hg 



Resistance coefficient for 0 diffusion 

A 

( 

min/ ml-mm Hg 


KS 

Solubility coefficient for oxygen 

mm Hg/ml 

RKC 

^T 

Rate constant for red cell destruction 

hr-1 

>|s* 

MCHC 

Mean corpuscular hemoglobin concentration 

gm/10 ml 



GUYTON 

NOTATION 

this 

REPORT 

DEFINITION 

UNITS 

*M02 

M02,. 

tis 

Rate of oxygen utLllzatLon of renal cells 

ml/min 

** 

P02 , 

art 

Oxygen tension in arterial blood 

mm Hg 

POT 


Oxygen tension in renal tissue 

mm Hg 

POV 

P02 

veti 

Oxygen tension in venous blood 

mm Hg 

VP 

. PV 

Plasma volume 

1 

** 

Q02 ^ 
art 

Oxygen flow rate in arterial blood 

ml/min 

*DOB 

Q02,. 

tis 

Oxygen delivery rate to renal tissue 

ml/ rain 

** 

Q02 

veil 

Oxygen flow rate in renal veins 

ml/min 

RC2 

ROD 

Red cell destruction rate 

l/min 

RCl 

RCP 

Red cell production rate 

i/min 

VRG 

RCV 

Red cell volume 

1 

OSA 

S02 ^ 
art 

Fractional saturation of ojgrgen in arterial blood 

ND 

+OSV 

S02 

ven 

Fractional saturation of oxygen in renal vein 

ND 


TOM 

Tissue PO effect on oxygen utilization 

ND 



GUYTON 

NOTATION 

this 

REPORT 

DEFINITION 

UNITS 

- 

V02""f 

art 

Maximum oxygen flow rate in arterial blood 

ml/ min 

*Q02 

V02,. 

tis 

Total oxygen volume dissolved in renal tissue fluid 

ml 

** 

V02,. 

tis 

Rate of oxygen accumulation in renal tissue fluid 

ml/ min 


* Guyton computes this variable for non-renal, non-muscle tissue and therefore it. is not equivalent 
to the corresponding variable in a renal model. 

** Not computed explicitly in Guyton model. 

No comparable term in Gu5rton model. 


ND- No dimension 



FIGURE 5 

COMPARISON OF RED BLOOD CELL CONTROL MODELS 


Process 

Myirea and Abbrecht 

Hodgeson 

Guyton 

Oxygenation 

P02<atm) 

^--^>P02(art) — >S02{art) 

input: S 02 (art) * 

Pulmonary Model 

1 ->S02(art) 

Blood Flow to Kidney 

Not a model parameter; 
Considered to be constant 

Input: Renal Blood Flow* 

Circulatory Model 

1 — ►non-renal blood flow 

Oxygen Extraction 

Not Considered 

Input: Metabolic Rate 

a) O 2 Diffusion * 

b) Variable Metabolic Rate * 

Stimulus for EP 
Release 

"Effective” [HbO] (art) 

= [HbO] X viscosity factor' 

Renal Tissue PO 2 * 

Non-renal Tissue PO 2 

Plasma [EP] 

T (clearance); Vd^p 

Mean life span,Tpp; Vdpp 

Not Considered * 

RBC Production Rate 

a) P = Pppjd “ + Pq’ 

b) Maturation time « f( EP ) 

,P «= f (log CEP]) 

^ “ ^gain^ ^^bias ' ^®2,tis^* 

RBC Destruction Rate 

RBC Life Span 

D « RC/ Mean life span 

D “ RBC ^ "^clearance * 

Distribution in PV 

Input variable: PV * 

Input variable: PV * 

Fluid Shift IVlodel ►PV 


* DENOTES PROPOSED USE IN PRESENT 
PO 2 = Oxygen Tension 
SOg = Oxyhemoglobin Saturation 
HbO = Oxyhemoglobin Concentration 
[EP] = Erythropoietin Concentration 


MODEL 

Vd^ = Volume of Distribution of EP 

PV = Plasma Volume 

T = Clearance Constant 

P == Production Rate 



















Oxygen Extraction by Tissue (2/3; 3/3&4; 4/2-14) 


Although recognizing the importance of oxygen extraction in determining the 
amount of oxygen available to the tissue, Mylrea and Abbrecht do not con- 
sider this factor, per se, in their model because of "inadequate data". How- 
ever, they fomd that their simulations did not agree with experimental 
responses to hypoxia until they added a "viscosity effect" to their model. The 
stated purpose of using this term was to account for a reduction in blood flow 
(and hence a decrease in oxygen delivery) due to increasing hematocrit. How- 
ever, since blood flow was not used explicitly in their model the actual effect 
was to obtain an "effective oxyhemoglobin concentration" in arterial blood 
lower than the true value. The overall result is a reduction in the amount of 
ojiygen being delivered to the tissues. This same effect was achieved more 
realistically in Hodgson's model by including a term for tissue oxygen consump- 
tion and in Giyton's model by accounting for both tissue metabolism and trans- 
capillary oxygen diffusion. The Guyton model allows for a variable rate of 
metabolism and diffusion resistance. 


Stimulus for Erythropoietin Release and Plasma |Epj (2/4&S; 3/5-8; 4/l5) 

Mylrea and Abbrecht postulate that their "effective oxyhemoglobin concentra- 
tion" is the stimulus for erythropoietin release. The relationship between 
these two variables is a monotonically decreasing function derived from experi- 
mental data on small animals. A time delay, also based on experimental 
observation, accounts for the delay between changes in (HbO) and the associ- 
ated erythropoietin release. Hodgson postulates that the renal tissue oxygen 
partial pressure is the stimulus for erythropoietin release. It is assumed 
that venous-capillary PO is identical to tissue PO and this term is computed 

A ^ 

by comparing rates of oxygen delivery to, and oxygen extraction by, the renal 
tissue. Guyton does not account for erythropoietin release, per se, but uses 
the tissue oxygen partial pressure as the stimulus for erythrocyte production. 

An inverse linear function relating tissue PO^ to er 3 dhrocyte production is assumed. 
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Erythrocyte Production in Bone Marrow (2/6; 3/9; 4/16) 


In the Mylrea and Abbrecht model, hemoglobin production rate (which is 
assumed to be similar to er 5 ^hrocyte production) is postulated to have a small 
basal value and asymptotically approaches a maximum value as erythropoietin 
concentration increases. The model incorporates a normal erythrocyte matur- 
ation time delay which decreases from 4 to 2. 5 days as erythropoietin concen- 
tration increases. Hodgson takes er 5 d;hrocyte production rate to be a linear 
function of In which is similar to the dose response curve proposed in 
Mylrea and Abbrecht 's model, but based on different data. Guyton assumes 
red cell production to be proportional to the difference between 02 {ygen demand 
and supply (as measured in terms of ojgrgen partial pressure). This "drive 
factor" for erythrocyte production has a minimum value and is multiplied by 
a gain factor which describes the overall sensitivity of the erythropoietin 
release - erythrocyte production circuit. In spite of the fact that Guyton does 
not account for an erythropoietic stimulating factor in his model, his overall 
schema of relating tissue ojsygen supply to red cell production is quite similar 
to those adopted by both Mylrea and Abbrecht and Hodgson. 

Red Cell Destruction and Distribution (2/7-12; 3/l0&ll; 4/17-21) 

It was assumed in Mylrea and Abbrecht's model that the red blood cells have 
a finite life span and that the cells destroyed at any given time are distributed 
in a normal manner about a mean value. Total circulating red blood cell 
volume was calculated from the net total circulating hemoglobin by using the 
appropriate hematologic indices. Hodgson accounted for red cell destruction 
by multiplying the production rate of red cells by the mean red cell life span 
to give the net number of red cells released. Red cell destruction, in the 
Guyton model, is accounted for by using a clearance time constant (based on 
the' mean life span) and multiplying this by the present red cell volume to give 
the rate of destruction. The net red cell volume is determined by integrating 
the difference between the rates of production and destruction. In all three 
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models the hematocrit of circulating blood is computed by dividing the present 
red cell volume by the total blood volume. Plasma volume is an input vari- 
able in Mylrea and Abbrecht and Hodgson's models while Gu 3 d:on computes 
plasma volume from an intercompartmental fluid shift model. In the red cell 
production block (Figure 4), plasma volume is considered to be an input term. 

It should be. apparent that these models are not exactly alike. For example, . 
Guyton does not explicitly consider an erythropoietic stimulating agent while 
Mylrea and Abbrecht and Hodgson use experimental data to compute release 
rates, clearance rates, and plasma concentrations of erythropoietin. On the 
other hand, Mylrea and Abbrecht consider renal blood flow to be constant while 
it appears that a variable flow rate can have important consequences on red 
cell production in the models of Hodgson and Guyton. 

The most important test of a model is its ability to predict experimental obser- 
vations using parameters that have physiological significance arid using values 
for those parameters that are biologically plausible. Validation was only 
attempted for one of these models, that of Mylrea and Abbreeht's, and that 
was for the case of mice under chronic hypoxic conditions . An example of 
the results they obtained are shown in Figure 6. It can be seen that the pre- 
diction of several important hematologic variables was in good ^reement with 
experimental data. Neither Hodgson nor Guyton have reported any attempt to 
subject their models to similar verification. 

New Model Proposed by Guyton-White 

Recently Drs. Guyton and White have been attempting to formulate an improved 
version of their erythropoiesis model. Their new model is quite preliminary 
and still untested, and it is described here briefly in the interest of complete- 
ness. It is based on a new concept of the renal tissue site that is believed to 
sense oxygen supply and regulate the release of erythropoietin. 
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Tlnne at 360 mmHg (days) 


Experimental data (Mylrea & Abbrecht, 1970) and the model results for mice 
exposed continuously to a total pressure of 360 mmHg (po 2 — 75 mmHg). Each point is 
the average of values obtained for nine mice. Different symbols represent different runs. 
The solid lines were obtained by first doing a linear interpolation to obtain daily values 
for the variables for each run, and then averaging the daily values for all runs at that 
altitude. Broken lines show the model output for this altitude, with the exception of the 
lowest curve, which shows the empirically derived relationship used for plasma volume in 
the model. Experimental values shown for time “0” were pbtained on control mice kept 
at atmospheric pressure. 


FIGURE 6 
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They have ascribed several characteristics to this oxygen sensing-erythro- 
poietin releasing tissue within the kidney: 1) the blood flow rate to this 
tissue is assumed constant due to excellent intrarenal autoregulation, 2) the 
oxygen consumption rate is also assumed to be constant, and 3) the erythro- 
poietin releasing system is considered to be a fast responding, high gain 
system, especially in hypoxic situations. The constancy of blood flow and 
tissue metabolism in effect make this tissue a sensor of blood hematocrit for 
the reasons that follow. Since oxygen delivery to tissue is dependent on basi- 
cally two parameters, the blood flow rate and hemoglobin concentration in the 
arterial blood (proportional to the hematocrit), a system that can keep blood 
flow constant will deliver ojgrgen at a rate proportional only to the hematocrit. 
In addition, by requiring a constant rate of ojQrgen consumption, the tissue PO 

U 

will vary in direct proportion to the hematocrit. It is postulated that erythro- 
poietin release is inversely proportional to hematocrit and is much more 
sensitive to charges in hematocrit below normal than above normal. 


There are some good arguments for an erythropoietic control system that is 
not dependent on blood flow. Primarily it is felt that, in most normal circum- 
stances', a change in tissue oxygen supply that is not in balance with demand 
can be compensated for by a rapid adjustment in blood flow and perhaps blood 
oxygenation by respiratory and circulatory reflexes. The more slowly respon- 
ding erythropoietic system is only required to be brought into play in those 
more miusual situations such as sustained breathing from hypoxic atmospheres 
that require additional correction. 

GE Model for Control of Erythropolesis 

The model represented in Figure 4 is being proposed as a starting point for 
simulating the human erythropoietic system. In particular we desire to inves- 
tigate the decrease in red cell mass that is observed in astronauts upon their 
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return from mission and the subsequent red cell mass increase following 
recovery. The model will be used to test several hypotheses, discussed in 
the next section, to explain this erythropoietic response. 

The proposed model is, in large measure, derived from a portion of the 
Guyton-White formulation (13, 14). Two segments of that model have been 
removed for study on a stand-alone basis: a) computation of net red cell 
mass from subroutine "BLOOD", and b) computation of tissue PO from 
subroutine "AUTRG". Some obvious changes were necessary in extracting 
a section of a larger model. Primarily these changes were concerned with 
input/ output flow. Other changes were made within the model itself in order 
to make it more physiologically plausible. These modifications are outlined 
below: 

1. Plasma volume (PV), blood flow (BF), arterial blood oxygen saturation 

(S02 ) and the autonomic stimulatory effect on metabolism (AOM) are 

computed from other physiological systems in the Guj^on model, but they 
will be considered to be input terms in the present model. 

2. The use of "non-muscle, non-renal" tissue oxygen tension as the primary 
stimulus for eiythropoiesis in the original Guyton model does not appear 
to be consistent with present physiological concepts. Consequently, the 
proposed model will assume that blood flow, tissue oxygen consumption, 
and tissue oxygen tension are variables associated with the kidney rather 
than with non-muscle, non-renal tissue. 

3. An algorithm has been included in the present model that computes per- 
cent saturation of-blood hemoglobin from the blood oxygen partial pressure 
and vice-versa. This algorithm is based on ah empirical formulation of 
the entire oxy-hemoglobin saturation curve (15). The Guyton model com- 
putes these terms by assuming a constant slope of the saturation curve in 
the appropriate range of oxygen tension. 
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Detailed Description and Assimiptions of the GE Model 


The flow chart of Figure 4 will be described more completely, especially 
those segments that differ from the Gu 3 d;on model (14). 

The computation of venous PO is described in Blocks 1 through 4. Block 1 

^ • 

xnd,x 

computes the maximum ojcjrgen flow possible, V02 , in the blood for a 

given hematocrit (HCT) and blood flow (BF). The mean corpuscular hemo- 
globin concentration (MCHC) and the oxygen carrying capacity of the hemo- 
globin molecule (CHbO } are considered constant. The fraction of the maximum 
oxygen flow that is delivered to the tissues (F02 ) is computed in Block 2. 

lIS 

The difference between the fraction of maximum' oxygen flow in arterial blood 
(S02 ) and the fraction of maximum oxygen flow diverted to the tissues 

CLI^h 

(Block 3) is equal to the fraction of maximum oxygen flow remaining in the 
venous blood Block 4 represents the oxyhemoglobin dissociation 

curve used to compute venous PO . 

The computation of tissue PO (Blocks 5 through 14) is quite similar to the 

development in the Guyton model. It is postulated that the "oxygen detector" 

site is located at the venous end of the capillary and therefore capillary PO 

is identical to venous PO This is derived from the assumption that the 

erythropoietin releasing tissue should be highly sensitive to oxygen lack and 

that small changes in blood oxj^en pressure would effect the venous PO^ to a 

relatively greater extent than arterial PO Tissue oxygen consumption 

(M02 ), Block 9, has a basal level, but can change due to: a) overall 

uxs 

autonomic stimulation (AOM), and b) an inhibitory effect caused by insufficient 
oxygen supply (TOM). 

The calculation of net red cell production rate is computed in Blocks 15 
through 17. This segment of the model is nearly identical to the Guyton for- 
mulation, Erythropoietin production, per se, is not considered in this 
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preliminary version of the model inasmuch as Information is not available 

regarding changes in plasma erythropoietin during space flight. The drive 

factor, DF (Block 15) represents the difference between tissue oxygen demand 

and supply in terms of oxygen tension. The gain factor {Block 16) is 

G 

merely the slope of a dose - response curve relating tissue PO to red cell 
production. Since the source of the data for generating the Guyton er 3 rthro- 
cyte production curve (Blocks 15 and 16) is not apparent, the proposed model 
will utilize a dose-response curve derived from the data of Hodgson (see Figure 
7). This is an extremely important element in the model because it represents 
the system that directly controls erythrocyte production. It should be recog- 
nized that the controlled variable is not red cell production or red cell mass, 
but a related term, tissue PO^. This is consistent with the accepted concept 
that the red-cell-hemoglobin mass is maintained at such a level that the 
tissue requirements for oxygen are optimally met. 

The calculation of hematocrit from the instantaneous red cell mass is des- 
cribed in Blocks 19 throi^h 21. Plasma volume is an input variable. This 
completes the analysis of the closed loop system of Figure 4. 


This model has recently been programmed for interactive -demand terminal 
use on the Univae 1110. An example of the response of the model, illustrating 
the associated graphics capability, is presented in Figure 8. 

4.0 CONCLUSIONS 

If the systems analysis presented here is basically correct (even though some 
of the details will be improved), it provides us with a systematic approach to 
investigate the heretofore unexplained changes in red cell mass observed in 
the Skylab program. For example, the analysis has identified five primary 
input variables, each of which could have a crucial effect on the response of 
the erythropoietic system: blood flow, plasma volume, mean red cell life 
span, sensitivity of the erythropoietin-erythrocyte producing circuit, and the 
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EFFECT OF TISSUE PO 2 ON RED CELL PRODUCTION 
(DERIVED FROM HODGSON, 1970) 

FIGURE 7 
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FIGURE 8 

Example of the computer generated graphics capabilities associated with the conversational mode of 
simulating the erythropoietic model of Figure 4. At t = 2 days (arrow) the arterial oxygen tension 
is changed from 95 mm Hg to 75 mm Hg, simulating a change in altitude from sea level to 5,000 ft. 
The sudden decrease in blood oxygen tension stimulates erythrocyte production which gradually de- 
creases to a new steady state value as the increased red cell mass is now able to provide sufficient 
tissue oxygenation. (Changes in oxygen-hemoglobin affinity, blood flow and plasma volume which 
are known to occur in hypobaric hypoxia were not considered in this example. 
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autonomic effect on cell metabolism. It is proposed to utilize digital com- 
puters to simulate the response of this model and to investigate the effects of 
the above mentioned parameters on the behavior of the circulating red cell 
mass. 

Hypothesis for Erythropoietic Adaptation to Weightlessness 

Based on this characterization of potentially important variables and taking 
additional physiological facts into consideration, it was possible to make some 
preliminary testable hypotheses as a starting point for investigating the de- 
creased red cell mass observed during space missions. It is postulated 
(Figure 9) that the stress of the zero-g environment decreases plasma volume 
and increases renal blood flow. These factors have already been implicated 
in the systems analysis as. potentiators for decreasing erythrocyte production. 

Preliminary data from the Skylab experiments strongly suggest that plasma 
volume decreases, at least early in the mission, leading to an increased 
hematocrit. It is certainly compelling to assume that this condition of "stress 
Induced hypovolemic polycythemia" results in a relative tissue hyperoxia and 
eventual shut-down of the erythrocyte producing circuit until the hyperoxic 
condition is relieved. But this has yet to be established in fact. Other un- 
known factors, induced by the still obscure physiological effects of weight- 
lessness, may be acting on the erythropoietin releasing-stem cell differen- 
tiation axis. 

The increase in renal blood flow during weightlessness has not been confirmed, 
but there is sufficient evidence to show that changing posture from an erect 
to supine position (which is often used to simulate the circulatory responses to 
weightlessness) can increase renal flow appreciably (11, 16). There are at 
least two reasons for such an increase; 1) cardiac output increases about 
20% in the supine position and the renal blood flow, which is normally about 
a fifth of the cardiac output, receives a share of this extra flow, and 
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FIGURE 9 

HYPOTHESIS TO EXPLAIN DECREASED RED CELL MASS DURING WEIGHTLESSNESS AND HYPEROXiA 

(STRESS & ADAPTATION STAGE) 






2) the normal centrally mediated reflexes which provide peripheral vaso- 
constriction and increased arteriolar resistance in the erect position are 
abolished to a large extent in the recumbent position and this results in as 
much as a 20% decline in renal vascular resistance.. {A mathematical model 
has recently been reported that is capable of simulating renal circulatory 
function during adaptation to changing gravitational forces (17) ). Further- 
more, it is known that this increased renal plasma flow during recumbency 
results in a less than proportional increase in glomerular filtration rate (16). 
This suggests that the recumbent renal flow is higher in all parts of the kidney 
vasculature including the vessels beyond the renal corpuscle. It is quite 
possible that under these circumstances the erythropoietin producing tissues 
are subject to an increased blood flow as well, leading to a relative hyperoxic 
state. Also of interest is the fact that renal flow is known to decrease as a 
result of exercise, even in. the recumbent position (18). This may help to 
explain the smaller decreases in red cell mass observed in those Skylab crew- 
members who were known to have engaged in the heaviest exercise regime. 

As suggested in Figure 9 an increase in atmospheric o^qrgen tension has been 
found to suppress the release of erythropoietin (19) as well as increase the 
rate of red cell destruction (1). While an increase in atmospheric oxygen 
tension was not a major factor in the Skylab missions, this hypothesis can 
explain the fact that in certain Gemini flights (an atmosphere of 100% oxygen 
at 258 mm Hg) the decline in red cell mass was greater than could be accounted 
for by a complete shutdown of erythrocyte production alone (2, 20). 

Recommendations for Future Study 

As a result of this review and analysis several areas of potential interest have 
been identified for more detailed study: 1) the response of renal blood flow 
to the stimulus of weightlessness, 2) the relationship of renal blood flow and 
production of erythropoietic factor, and 3) the relationship between the response 
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of the cardio-respiratory system and the ej^rthropoietic system to abnormal 
tissue oxygen balance. Where possible, the present literature survey will be 
expanded to include coverage of these items. 

In addition, current periodicals will be reviewed on a continual basis to keep 
abreast of new developments in the field of erythropoietic control especially 
as concerns identification of other chemical agents or neural components that 
affect red cell production and the clarification of the role of tissue oxygen as 
a stimulus for erythropoietin production. 

A model claiming to represent the control of erythropoiesis should be able 
to simulate not only the data obtained during the Skylab program, but should 
also be capable of simulating the responses to a wide range of hemotologic 
disturbances such as hypoxia, hyperoxia, hemorrhage, anemia, polycythemia, 
and hemolytic diseases, as well as responses to immersion and bedrest. For 
this reason, a literature survey will be made to locate data regarding these events. 

Many factors have been discussed that are known to affect erythropoiesis, but 
have not been included in the proposed model. Performing computer simula- 
tions of a mathematical model is not unlike performing experiments in the 
laboratory on an unknown system, starting from asking simple and easily 
imderstood questions and progressing to the more complex world of reality. 

It is probably a mistake to try to start from too realistic a system. Accord- 
ingly, it is believed that at least some, if not all, of the following factors may 
eventually be necessary in a more complete model: 1) a description of an 
erythropoietic stimulating factor (ESF), 2) a description of a still-not-con- 
firmed second erythropoietic factor that is perhaps a precursor to ESF, 

3) a description of the kinetics of stem cell differentiation in the bone marrow, 

4) a more accurate description of the renal sensor slte(s) of oxygen tension, 

5) consideration of a possible shift in the affinity of hemoglobin for oxygen, 
and 6) consideration of possible extrarenal sites of ESF production. 
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Assessment of the Systems Analysis Approach 

It is believed that the utility of the systems analysis approach as an investi- 
gative tool in the Sl^lab hematological experiments may be limited by the 
scarcity of measurements of certain potentially crucial parameters, as well 
as an incomplete understanding of the mechanisms affecting erythropoiesis. 
However, it is felt that mathematical modeling and simulation can be quite 
useful in several ways; 1) identification of important parameters and their 
sensitivity on the overall system, 2) a method to test hypotheses rapidly, 

3) identification of specific elements that must be further quantified, and 4) 
prediction of the behavior of certain unmeasured or difficult to measure para- 
meters of the system in response to specific stress states. 

There is every reason to suspect that this study, in its initial stages, will 
suggest more questions than answers. However, this should help point the 
way toward the kind of clarifying experiments that might be performed and the 
type of data that should be collected. As the simulation study progresses it 
may be possible to incorporate more and more diverse kinds of experimental 
results into a single model, each of which alone would not support a generalized 
theory, but taken together might point to a coherent picture of erythropoietic 
control. 
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